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Abstract 

The branching ratio of beauty hadrons to final states containing two charm hadrons, 
Br(b DDX), has been measured using an inclusive method in hadronic decays with the 
OPAL detector at LEP. The impact parameter significance of tracks opposite tagged b-jets is 
used to differentiate b — > DDX decays from other decays. The result is 

Br(b ^ DDX) = (10.0 ± 3.2(stat.)l2:9(det.)+^V(phys-))%> 

where "det." is the systematic uncertainty due to the modelling of the detector, and "phys." is 
the systematic uncertainty due to the modelling of the underlying physics. Using this result, 
the average number of charm plus anti-charm quarks produced in a beauty quark decay, ric, is 
found to be l.UtoH- 
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1 Introduction 



Studying the decays of b hadrons allows important tests of the Standard Model and Heavy 
Quark Effective Theory (HQET) to be made. One such test is whether the average number of c 
plus c quarks produced in the decays of b quarks, ric, is consistent with theory. The quantity ric 
can be determined experimentally by measuring the "topological" branching ratios of b hadrons 
to different numbers of charm hadrons: 

nc = 1 + Br (b DDX) + Br(b — >^ charmonium) — Br(b — >^ no charm). (1) 

This analysis measures the inclusive branching ratio of b hadrons to two charm hadrons, 
Br(b— *^DDX), and combines this number with previous measurements of Br(b — > charmo- 
nium) and Br(b no charm) to obtain ric- Neubert and Sachrajda use HQET to calculate 
ric = 1-20 lb 0.06 [1]. This theoretical prediction for Uc is currently limited by uncertainty in the 
ratio of the charm and beauty quark masses (0.25 < mc/m;, < 0.33). 

Besides being interesting in their own right, Br(b — *■ DDX) and ric are correlated to the b 
hadron semi-leptonic branching ratio, Br(b ivX.). The current combined experimental values 
for Br(b ^ luX) = (10.59 ± 0.22)% at the Z° (LEP) [2] and Br{B^/+ luX) = (10.38 ± 0.32)% 
at the T(4S) (CLEO) [2] are slightly lower than the central values of theoretical predictions. 
QCD calculations within the parton model yield Br(b — > £i/X) > 12.5% [3,4], while more recent 
calculations that include radiative QCD corrections, spectator quark effects and charm quark 
mass effects yield Br(b ii'X) = 9.5% 13.0%, depending on the renormalization scale and 
the quark mass scheme (pole or MS) used for the calculations [1,5]. If any component of the 
hadronic width {e.g. h DDX) is larger than expected, then the central value for Br(b £vX) 
will be smaller. 

The analysis presented in this paper makes the first inclusive measurement of Br(b DDX) 
and ric using OPAL data. This analysis uses a technique similar to one employed by DELPHI [6] : 
a joint probability variable, constructed from track impact parameters, is used to discriminate 
amongst the different b hadron decay topologies. 

2 The OPAL Detector and Data Samples 

A brief description of the most relevant components of the OPAL detector is given here; see 
Reference [7] for more details. The central tracking system consisted of a silicon microvertex 

detector, a precision vertex drift chamber, a large volume drift (jet) chamber, and a set of 
chambers surrounding the jet chamber that made precise measurements of the z-coordinates of 
tracks^ . The silicon microvertex detector consisted of two layers of silicon strip detectors that 
provided three-dimensional hit information from 1993 onward. The tracking system was located 
inside an axial 0.435 T magnetic field generated by a solenoidal coil just outside the tracking 
chambers. Surrounding the solenoid were, in order, the scintillation time-of-flight detectors, 
the pre-sampling devices for the electromagnetic calorimeter, the lead glass electromagnetic 
calorimeter, the iron return yoke for the magnetic field (instrumented in order to provide hadron 
calorimctry), and farthest from the interaction point, the muon chambers. 

The data used in this analysis were collected from 1993 to 1995. During these years, the 
centre-of-mass energy of the colliding beams, y/s, was approximately mzo. Data collected off 

'^In this paper, Roman font "b" refers to the admixture of weakly decaying hadrons containing a beauty quark 
produced in e+e~ annihilations at ^/s = in^^o; italic font "6" refers to the beauty quark. Weakly decaying hadrons 
containing a charm quark, c, are collectively called D hadrons. 

^The OPAL coordinate system is right handed with the «-axis following the electron beam direction, and the 

X-axis pointing to the middle of the LEP ring. The azimuthal angle <f> is measured in the x — y plane with respect 
to the +x-axis and the polar angle 9 is measured with respect to the +z-8ixis. 
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the Z peak are not used because the joint probabihty distributions used in this analysis vary as 
a function of centre-of-mass energy. After applying standard OPAL hadronic TF^ decay selection 
cuts [8], 1,866,000 events remain for analysis, with the contribution from background less than 
0.1%. 

A total of 8 million simulated hadronic li^ decays are used to generate probability density 
functions (PDFs) for all signals and backgrounds considered for this analysis. Of these decays, 
3 million are Z*^ bb decays. The simulated hadronic decays are generated by JETSET 7.4 [9], 
using parameters tuned by OPAL [10], then processed using a full Monte Carlo (MC) simu- 
lation of the OPAL detector [11]. Both real and simulated events are subjected to the same 
reconstruction and analysis algorithms. 



3 Method 

3.1 Joint probability 

The main goal of this analysis is to differentiate double charm b hadron decays from single 
charm b hadron decays. Single and double charm decays can be statistically separated by their 
different topologies. In b-jets, tracks from D hadron decays originate farther from the interaction 
point (IP) than tracks from b hadron decays. Most of the tracks in b ^ DDX decays originate 
from D hadron decays, while most of the tracks in single charm b hadron decays originate from 
the location where the b hadron decays (see Figure 1). The separation between a track and the 
IP (approximated by the reconstructed primary vertex) can be expressed in terms of its signed 
I — (f) impact parameter significance^, S = do /(J do i where do is the impact parameter of the track 
with respect to the primary vertex in the r — (/> plane, and (T^q is the uncertainty in this quantity. 
The S values of the tracks in a jet are used to calculate a single variable for each jet: the joint 
probability, Pj. 

The joint probability is calculated by first considering the S value of each selected track in a 
jet. Under the hypothesis that each track originates at the IP, one can calculate the conditional 
probability, p^, for a track with 5 > to have the measured value S or larger. This probability 
is calculated by integrating the measured S resolution function for the OPAL detector, f{S), 
beyond the measured value of S: 

f^^"* fiS)dS 

fis)ds 

where S'cut (= 25) is a cutoff in S beyond which tracks are not considered. Given an ensemble 
of tracks that originate from the interaction point, the distribution oi pi will be uniform from 
to 1. 

The joint probability, Pj, is calculated by considering the pi of all the tracks in a jet: 

^-^ ml 

where y is the product of the N individual track probabilities. Pj is the probability that the 
product of A'^ random numbers uniformly distributed from to 1 is y or smaller. The larger the 
S of tracks in a jet, the smaller the pi, and hence also Pj, will be. The possible values of Pj go 
from to 1, so — In(Pj), the variable used in the analysis, varies from to +00. 

^The sign of S is determined by the location in the r — (j> plane where the track crosses the jet axis. The sign is 
positive (negative) if a line drawn from the primary vertex to the point where the track and jet axis cross points 
in the same (opposite) direction as the jet. 
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b decay 

IP • 



charm 



IP • 



IP • 




1 charm 



2 charm 



Figure 1: Three different topologies of b hadron decays. As the charm hadron multiphcity 
increases, the tracks tend to originate farther from the interaction point (IP). Both b ^ no charm 
and b — > charmonium decays are represented by the upper "0 charm" diagram. These two 
decay modes have similar topologies because of the prompt electromagnetic or strong decays of 
charmonium states. 
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Br(b — > DDX) is measured by comparing simulated — In(Pj) distributions of different b 
hadron decay topologies and backgrounds to data. The data distributions of — In(Pj) are fit 
with the simulated — In(Pj) distributions with Br(b DDX) as a free parameter. 

3.2 Event, jet and track selection 

A number of event level cuts are applied to the data. The OPAL hadronic decay selection is 
applied and the silicon microvertex detector is required to have been fully operational. Events 
are also required to be well contained in the central portion of the tracking detectors (especially 
the silicon microvertex detector) by imposing a cut on the direction of the event thrust axis: 
I cos0t\ < 0.7. In addition, events are required to be two-jet-like by applying a cut on the thrust 
of the event: T > 0.85. The thrust and thrust axis of each event are calculated using tracks and 
energy clusters in the electromagnetic calorimeter that are not associated with any track. 

All events are divided into two jets by the Durham jet finding algorithm [12]. The OPAL 
LEP-2 b-tagger [13] is applied to these jets to select a sample of jets that is enriched in b 
hadron decays. This b-tagger has a good b-tagging efficiency while maintaining a high b hadron 
purity; this is achieved by combining track, high transverse momentum lepton, and jet shape 
information into a single likelihood variable. For this analysis, the b-likelihood is required to 
be greater than 0.9 in order to obtain a high purity b-jct sample (purity ~95% and efficiency 
~40%). A jet is used in the analysis if the opposite jet passes the b-tag cut. If both jets pass 
the b-tag cut then both jets are used. Using the opposite jet for b-tagging provides an unbiased 
sample of b-jets. 

Track selection cuts are made to select a sample of well measured tracks enriched in decay 
products of b and D hadrons. As a preliminary track selection, tracks are required to have an 
T — (f) impact parameter with respect to the beam spot of less than 5 cm, at least 20 jet chamber 
hits, momentum p < 65.0 GeV/c, and transverse momentum px > 0.15 GeV/c. In addition, 
tracks are required to have r — and z coordinate hits in both layers of the silicon microvertex 
detector as MC studies show this requirement greatly reduces the systematic uncertainty due 
to detector resolution modelling. In order to reduce the number of fragmentation tracks, tracks 
are required to have a signed impact parameter significance, 5 > 0, an angle with respect to the 
jet axis, dt-j < 0.6 radians, and a rapidity with respect to the jet axis, y = \ ln( g^^|| ) > 1.0 {E 
is the energy of the track and p|| is the component of the track's momentum parallel to the jet 
axis). The efficacy of a p cut for reducing the number of fragmentation tracks was investigated; 
however, MC studies show that a p cut (beyond the px > 0.15 GeV/c requirement) resulted 
in an unnacccptably large increase in the statistical uncertainty of Br(b DDX). The overall 
selection efficiency for tracks from b and D hadron decays is approximately 40%. 

3.3 S resolution function determination 

In order to calculate Pj, the S resolution function of the OPAL detector, f{S), must be known. 
The S resolution functions are determined by fitting functions that are the sums of three gaus- 
sians plus an exponential, to tracks with 5 < (backward tracks) in a heavy flavour suppressed 
data sample (opposite jet has b-likelihood < 0.1). Backward tracks are used as they tend not to 
have genuinely positive impact parameters; the S distribution of backward tracks is dominated 
by detector resolution effects. Requiring b-likelihood < 0.1 reduces the fraction of selected back- 
wards tracks from b or D hadron decays from 11% to 5% (4% from D decays). The presence of 
a small fraction of tracks from b and D hadron decays is not critical because all that is required 
is a sample of tracks to determine a common function, f{S), that will be used to calculate pi 
and Pj for data and MC. The same event and track selection cuts are applied to select backward 
tracks for f{S) and "forward" tracks for Pj (except by definition, S < for backward tracks 
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and 5' > for forward tracks). 

For each year of data taking, different f{S) are determined for tracks in three momentum 
bins: p < 1.5 GeV/c, 1.5 GeV/c < p < 4.0 GeV/c, and p > 4.0 GeV/c. The agreement between 
simulated and real data for the S distributions of backward tracks is shown in Figure 2. The 
agreement is good for S values less than 5 (where 99.7% of the tracks lie); however, there 
are some differences in the high S value tails of the distributions. The systematic uncertainty 
associated with this disagreement is discussed in Section 5.1. 

3.4 Fitting procedure 

A fit is performed to estimate the best fit values of Br(b — > DDX) for each year's data. The 
MINUIT package [14] is used to minimize the function. The fitting function used is 

F{x) = N{1+ ax) 



{1 fvds fc fg} ^ 

{BrocG°'(x) + BricG^'^(x) + Bi2cG^^{x) + Bi^G^{x)} + 
{fudsG^^'^x) + hG\x) + fgG^{x)] 



(4) 



where x = — In(Pj) and AT is a normalization factor chosen so that / F{x)dx is equal to the 
number of events in the data. The G^{x) are the normalized PDFs for the different signals and 
backgrounds. Bric is the b single charm branching ratio; Br2c is the double charm branching 
ratio; Broc is the no charm branching ratio; Br^ is the b^ charmonium branching ratio. The 
background fractions are fc, the fraction of charm jets, fuds^ the fraction of light quark jets, and 
fg, the fraction of gluon jets from Z° bb events. Finally, a is a term used to parameterize 
mis-modelling due to incomplete knowledge of all physics inputs. 

MC studies show that changing the input values for various physics inputs such as the mean 
multiplicity of charged particles from fragmentation ((A'ch)frag)) b hadron lifetimes (r^) or a 
variety of other inputs results in approximately linear changes in the — In(Pj) distributions. The 
inclusion of the a term in Equation 4 reduces the sensitivity of the analysis to changes in these 
physics inputs as this term allows the data to constrain these inputs. Consequently, several of 
the main systematic uncertainties are reduced by up to 50% compared to using a fitting function 
without a. This results in the total uncertainty being reduced by 30%. 

Only Br2c and a are free parameters in the fit. Bric is given by 

Bric = 1 - Br2c - Broc - Br^,, (5) 

while Broc = (0.7 ± 2.1)% [6], and Br^ = (2.4 ± 0.3)% [15] are fixed. The background fractions 
fuds, fc and fg arc also fixed in the fit; the values of the fractions are determined from the MC. 

The value of a is constrained to be close to zero by including a "penalty function" in the 
evaluation of the so that 



a 



2 



X'->X' + — • (6) 

MC studies were performed to determine the value of cTq in Equation 6 that results in the 
smallest predicted total uncertainty for Br(b DDX): (Tq = 3 x 10~^. This magnitude of 
o"a corresponds approximately to the change in a when the sources of the main systematic 
uncertainties arc varied by their one standard deviation uncertainties. The relative sizes of 
the statistical and systematic uncertainties changed appreciably as cr^ was varied in the MC 
studies; however, due to the anti-correlation between these uncertainties, the total uncertainty 
for Br(b — DDX) was not very sensitive to the value of aa- Note that MC studies show that 
introducing the a parameter in the fit does not bias the estimator for Br(b DDX) and reliable 
statistical uncertainties are estimated. 
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Figure 2: Impact parameter significance distributions for backward tracks with hits in both 
layers of the sihcon microvertex detector. The plots in the left column show backwards tracks 
with < —S < 5. The data points with error bars represent the 1994 data and the histograms 
represent the corresponding simulation. The plots in the right column show the ratio of MC/data 
for S values < — < 15. 
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3.5 Binning of data by track multiplicity 



The data are divided into six track multiplicity (TM) bins to improve the sensitivity of the anal- 
ysis and to reduce uncertainties due to incorrect modelling of the number of tracks contributing 
to — In(Pj). There is one bin for each value of track multiplicity in the range 1 < TM < 5 and 
one bin for TM > 6. 

In order to ensure that the minimum number of expected data entries in each — In(Pj) bin is 
at least 50, the number and range of — In(Pj) bins varies for each track multiplicity. Jets with 
—\n{Pj) greater than the range considered are included in the largest —\n{Pj) bin. The first five 
— In(Pj) bins of each track multiplicity bin (ten bins in the case of TM=1) are also combined into 
one large bin in order to reduce the sensitivity of the analysis to changes in the do resolution. 
Figure 3 shows the MC —In(Pj) distributions for single and double charm b hadron decays for 
the different track multiplicity bins. 

Because the data are binned by track multiplicity, it is the fraction of b — > DDX decays in 
b hadron decays for a particular track multiplicity bin, /2cj, that is actually determined in the 
fit for each bin. According to the MC, f2ci increases by a few percent (absolute) as the track 
multiplicity increases from one to six. Br(b DDX) is determined by summing the results for 
all track multiplicity bins. Using the fitted values of f2ci for each bin and N^. , the total number 
of b hadron decays in each bin, the number of b ^ DDX decays (= f2c. x Nh.) is determined 
for each bin. Br(b DDX) is calculated for each year's data by dividing the total number of 
b DDX decays by the total number of b hadron decays: 

Br(b ^ DDX) = ^^=1 ^'^^^^''i . (7) 
Note that a single value of a is determined for all track multiplicity bins in a year. 



4 Results 

4.1 Results for each year 

The results of the fits for each year are shown in Table 1. The probabilities to obtain the 
values in Table 1 or larger for 84 degrees of freedom arc 0.36, 0.16, and 0.59 for 1993, 1994 and 
1995 respectively. The fitted values of Br(b — > DDX) for each year and each track multiplicity 
bin are not constrained to the physically allowed region, 0% < Br(b DDX) < 100%, because 
the results for each year and each track multiplicity bin are combined [16]. 



year 


data events 


Br(b ^ DDX) (%) 


a (xlO-3) 


xVd.o.f. 


1993 


408k 


-2.2 ±6.5 


1.0 ±2.4 


88.0/84 


1994 


1,076k 


15.0 ±4.4 


-0.8 ±2.0 


96.6/84 


1995 


382k 


15.5 ±6.7 


-1.4 ±2.7 


80.4/84 



Table 1: Results of — In(Pj) fits for each year of data-taking. The uncertainties are statistical 
only. 



The best fit values of a are consistent with zero for each year of data taking. This shows 
that the physics inputs to the MC and the detector modelling are in reasonable agreement with 
the data. Due to differences in detector performance from year-to-year, it is not expected that 
the value of a should be the exactly same for each year; it is for this reason that a is determined 
separately for each year. Note that the central value of Br(b DDX) changes by less than one 
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Figure 3: The 1994 MC —ln{Pj) distributions for single and double charm b hadron decays 
for each track multiplicity bin. The last bin of each histogram includes all data beyond its 
upper edge. The solid and dashed histograms are for single and double charm b hadron decays, 
respectively. 
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standard deviation when a is either unconstrained or omitted from the fitting function in fits to 
data or to a distinct sample of MC used instead of data. 

Figure 4 shows the sum of the fitted simulated PDFs and data for all years. For TM>6, the 
single charm component is fitted to be > 100% and Br(b DDX) < 0%; for the correspond- 
ing plot, there is no visible b — > DDX component and the number of entries in the light grey 
histogram is greater than the number of entries in the line histogram (sum of all the MC PDFs). 

A cross check of the analysis is performed by repeating the analysis for 1994 on a MC sample 
instead of data. The Br(b DDX) result is (17.4l^'g(stat.))%, which is consistent with the 
true Br(b DDX) value in MC, 13.3%. The fitted value of a is (-0.02 ± 2.01) x IQ-^. 

4.2 Combination of Br(b DDX) results 

The Br(b DDX) results for each year of data taking are combined to yield 

Br(b ^ DDX) = (10.0 ± 3.2(stat.)l|^(det.))%, (8) 

where "det." is the uncertainty due to detector modelling. The uncertainty due to detector 
modelling is considered to be uncorrelated from year-to-year as the detector modelling was 
tuned separately for each year. The size of this systematic uncertainty is approximately the 
same for each year. Systematic uncertainties from the modelling of particle physics processes 
(sec Section 5.2) are fully correlated from year-to-year so do not need to be considered when the 
weighted mean for Br(b DDX) is calculated. The weights for each year are l/(o"stat. + ^dct )• 
The x^/d.o.f. for combining the results from the three years is 3.7/2. The probability to obtain 
this ov larger with two degrees of freedom is 0.16. The Br(b DDX) results can also be 
compared by track multiplicity bin. Note however, that it is expected that /2cj should be 
slightly different for each bin. The x^/d.o.f. for combining the results of each bin is 9.8/5. The 
probability to obtain this or larger with five degrees of freedom is 0.08. 

5 Systematic Uncertainties 

The sources of systematic uncertainty investigated for this analysis are discussed in Sections 
5.1 (detector modelling) and 5.2 (particle physics modelling). Tables 2 and 3 summarize the 
systematic uncertainties. 

5.1 Detector modelling 

Applying the same set of f{S) to both data and MC to calculate —\n{Pj) assumes that the S 
resolutions are the same between the two. The MC has been tuned to make the do and 
distributions of backwards tracks in the MC match the data distributions (see Figure 2). The 
MC do distributions are tuned by scaling the difference between true and reconstructed track 
parameters. The MC a^g distributions are tuned by scaling the values of cr^o by a constant. 

The uncertainty due to mis-modelling do in the MC is determined by re-processing the MC 
with the do resolution varied by ±1.5%. This variation accounts for differences between the 
tuned MC track parameters determined by three different tuning methods. The MC is tuned by 
comparing MC and data backwards track do distributions at small do (width of core gaussian 
describing do distribution) and over a wide range of do, and by comparing — In(Pj) distributions 
for backwards tracks in a jet. 

The uncertainty due to mis-modelling a^^^ is determined by repeating the analysis with a^^ 
for each track scaled by ±0.5%. This variation accounts for remaining disagreement between 
the data and the MC. 
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Figure 4: Sum of fitted MC PDFs and data for all years for each track multiplicity bin. The 
data points represent the data and the line histograms are the sums of the MC PDFs. The light 
grey histograms depict the single charm b hadron decays, and the hatched histograms depict 
the b — > DDX component. The dark grey histograms depict all backgrounds. The b DDX 
contribution for each plot is calculated from the fitted Br(b — > DDX) values for each track 
multiplicity bin from each year. 
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Sign of 


Source 


Value 


^Br(b^DDX) C^) 


ABr(b^DDX) 
ASource 


do modelling 


±1.5% 


+2.2 
-2.8 


N/A 


CTrfg modelling 


±0.5% 


±0.5 


N/A 


^track 


see text 


±0.2 


N/A 


Uncorrelated total 




+2.4 
-2.9 




(-^ch)frag 


12.46 ± 0.32 


±6.2 


+ 




see text 


±0.5 


N/A 


TbO 


(1.542 ±0.016)ps 


±0.9 




7"B+ 


(1.674 ±0.018)ps 


±1.0 






(1.461 ±0.057)ps 


+0.7 

-0.6 






(1.208 ±0.051)ps 


+1.0 
-0.9 




td+ 


(1.051 ±0.013)ps 


±0.1 






(0.412 ± 0.003)ps 


±0.3 






(0.490 ± 0.009)ps 


±0.3 






(0.200 ± 0.005)ps 


±0.2 






(10.5 ±2.0)% 


±1.7 


+ 




(9.2 ± 2.4)% 


+0.4 
-0.3 




g^bb 


(2.54±0.50)xl0-3 


±0.0 


N/A 


g cc 


(2.99±0.39)xl0-2 


±0.2 


+ 



Table 2: Summary of systematic errors for Br(b — > DDX) (part 1). The continuation of the 
summary, including the total correlated uncertainty, is found in Table 3. The definitions and 
explanations of all the sources of uncertainty are contained in Sections 5.1 and 5.2. 
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Sign of 


Source 


Value 


'^Br(b^DDX) ^/oj 


ABr(b-+DDX) 
ASource 


\"'cn/D+ 




+0.6 
-0.8 




(^ch)DO 




+0.9 
-1.3 




\"'ch/D+ 




+1.3 
-1.0 




\"'ch/\+ 


2.7 ± 0.5 


+1.1 
-0.9 




(t) n >T-v-L 
\'*'7r^/D+ 


1 18 ± 33 


+0.7 
-0.9 


-1- 




1 31 lb 27 


+6.4 
-3.5 


+ 
1 


\'^7rO/D+ 


9 0+14 


+1.8 
-1.4 






/fil 9-1-7 


+1.1 
-1.0 




XJi I -L' ^ -l-v -iV J 


MIX -U QW 


±1 3 




Br('D+ ^ K°X) 


(■39+28)% 


±1.9 






/ocj. 11 W 
\^0(J ^ ± ±^ /o 


+0 4 


_1_ 


\"'Cn/ D 


4 97 ± 07 


±0.7 




/d+(1c) 


(23.3 ± 2.9)% 


+1.3 
-1.4 




/d+(2c) 


(17.0 ± 4.9)% 


+0.3 
-0.2 




/a+(M 


(10.0 ± 2.9)% 


+2.3 
-2.4 


+ 


/a+(2c) 


(7.4 ±2.9)% 


+1.7 
-1.3 


+ 


^9 


±10% 


±0.4 


+ 




±10% 


±2.0 


+ 




±10% 


±0.8 


+ 


Br(b— >^ no charm) 


(0.7 ±2.1)% 


+3.0 
-1.0 


+ 


Br(b— ^ charmonium) 


(2.4 ± 0.3)% 


±0.3 


+ 


Correlated total 




+10.4 
-9.0 





Table 3: Summary of systematic errors for Br(b DDX) (part 2). The correlated systematic 
uncertainties are due to particle physics modelhng. 
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After tuning the MC do and distributions, a difference remains between the tails of the 
data and MC S distributions for backwards tracks (see Figure 2). A study of the impact of 
this difference on the measured value of Br(b DDX) was performed by varying, by 25%, the 
fraction of tracks in the MC whose measured do values are significantly different (greater than five 
standard deviations) from their true do values. The measured value of Br(b DDX) changes 
by a small amount (1.4%) that is already covered by the systematic uncerainties attributed to 
do and cr^,, modelling. 

Different track selection efficiencies, etracki in the MC and the data may result in systematic 
differences between the MC and the data joint probability distributions. The most significant cut 
for etrack is the cut on the number of silicon microvertex detector hits associated with a track. 
After correcting the fraction of tracks in the MC with associated silicon hits (by randomly 
dropping tracks with silicon microvertex detector hits from the calculation of — In(Pj)), only a 
small statistical uncertainty remains. 

5.2 Particle physics modelling 

The dominant sources of systematic uncertainty in this analysis are those which significantly 
affect the S values of tracks included in the calculation of ^ln(Pj). The sources of the largest 
uncertainties for Br(b DDX) arc the charged particle multiplicity from fragmentation in — >■ 
bb events, the neutral particle multiplicity of D decays, and the fractions of different D species 
in b hadron decays. Every correlated systematic uncertainty is calculated separately for each 
year, then combined for all years to yield the total systematic uncertainty due to each source. 
The total systematic uncertainty on Br(b DDX) due to physics modelling is calculated by a 
quadrature sum of the individually combined systematic uncertainties. 

Care must be taken when considering the systematic uncertainty due to fragmentation in 

— > bb events, as (xe), the mean energy of weakly decaying b hadrons in decays, is closely 
related to (A'ch)frag5 the average number of charged particles produced in the fragmentation 
process (i.e. charged particles not from the decay of the b hadrons). Reweighting jets to vary 
(xe) changes (A^ch)frag at the same time. The contribution of each needs to be separated in 
order to avoid double counting systematic uncertainties. 

In order to determine the systematic uncertainty due to (A''ch)frag independently of (xe), 
(A^ch)frag is Varied by randomly dropping fragmentation tracks. The mean multiplicity of charged 
particles from fragmentation in Z° — > bb events is determined by comparing experimental values 
of the average charged particle multiplicity in Z° bb decays, {nch)bi [24,25], and the average 
charged particle multiplicity of b hadron decays, (nch)b (including the charged decay products 
of Kg and A) [2]. Combining these measurements gives (iVch)frag = 12.46 ± 0.32. The difference 
between Br(b DDX) before and after dropping fragmentation tracks represents the systematic 
uncertainty. The process of randomly dropping tracks is repeated 20 times to obtain a more 
precise estimate of the systematic uncertainty. 

Several experiments have made precise measurements of (xe) [17-19]. The combined result 
for the model-independent value of (xe) is 0.7151 it 0.0025 [20]. To assess the uncertainty in 
Br(b DDX) due to uncertainty from (xe), jets originating from b hadrons in the simulation are 
reweighted so that the model-independent value of (xe) is varied by its one standard deviation 
uncertainty assuming the Bowler [21], Lund [22], and Kartvelishvili [23] fragmentation models. 

The change in Br(b DDX) due to changing (xe) was found to be consistent with being 
due entirely to the associated change in {Nch)fTa.g- This can be understood by realizing that 
the do of tracks from b and D hadrons do not change significantly as a function of (xe); as 
(xe) increases, the decay lengths of the b and D hadrons increase but the opening angles of the 
tracks with respect to the jets decrease. These two effects tend to cancel out. For this reason, 
no systematic uncertainty is attributed to (xe) itself. 
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CLEO has made precise measurements of the scaled^ momentum spectra of D'^*-' hadrons in 
B meson decays, (ajb-^o) [26]. Those results are applied to the admixture of b hadrons produced 
in decays. The uncertainty in Br(b — > DDX) due to the uncertainty of (xb^o) is determined 
by repeating the analysis many times with different, randomly generated MC (xb^o) spectra 
that are compatible with the CLEO data. The width of the distribution of Br(b — > DDX) 
values obtained with the different D momentum spectra determines the associated systematic 
uncertainty. 

The lifetimes of the b and D hadrons are independently varied by their experimental uncer- 
tainties quoted in [2]. The fractions of different b hadrons are also varied as prescribed by the 
LEP Heavy Flavour Working Group [27]. 

The g ^ bb and g ^ cc rates arc varied by their experimental uncertainties [27]. Requiring 
the thrust of events to be greater than 0.85 greatly reduces this background, so the resulting 
systematic uncertainty for Br(b — > DDX) is very small. 

The Mark III collaboration has published values for the mean number of charged particles, 
(nch)D, and neutral pions, (n^o)D, produced in D+, D° and D+ decays [28]. Br(D— K X) was 
also measured. These values, along with Br(A^ AX) [2] and the charged particle multiplicity 
in charm baryon decays (varied ±0.5 about the JETSET prediction), are separately varied in 
the MC by their uncertainties to determine their contributions to the systematic uncertainty 
of Br(b DDX). While each multiplicity is varied, the others are kept constant. Varying the 
neutral particle multiplicities affects the pT distribution of charged particles from D decays, 
which affects the S distribution of tracks contributing to — In(Pj). 

The dependence of Br(b DDX) on the mean charged particle multiplicity of b hadron de- 
cays, (nch)b) is reduced by binning the — In(Pj) distributions by track multiplicity, but there still 
exists a systematic uncertainty due to (nch)b [27]. The systematic uncertainty of Br(b DDX) 
due to the uncertainty of (nch)b is determined by reweighting jets to effectively change (rach)b 
in the MC. 

The fractions of different D hadrons in single charm, /^^(Ic), and double charm, foii^c), b 
decays are varied because the different D hadron species have quite different lifetimes: the ratio 
of Tj^+ : Tq+ : Too : is approximately 2.5 : 1.2 : 1 : 0.5. As the D"*" and A+ possess the 
longest and shortest D lifetimes, they are the two D hadrons considered in this section. Using 
measured rates of D hadron production in b hadron decays at the Z'^ [2] , the following fractions 
arc calculated: /d+(1c) = (23.3 ±2.9)%, /d+(2c) = (17.0 ±4.9)%, 4+(lc) = (10.0 ±2.9)%, and 
4+ (2c) = (7.4 ±2.9)%. 

The backgrounds are divided into the following categories: gluon jets in — > bb events 
(fg), light quark background (fuds), charm quark background (/c), b — no charm decays, and 
b charmonium decays. The fractions of backgrounds assumed to be present in the data are 
varied by their uncertainties to determine their contributions to the systematic uncertainty of 
Br(b — > DDX). The c and uds selection efficiencies in b-tagging are varied by ±10% [29]. The 
assumed branching ratios for b charmonium and b — ^ no charm are also separately varied to 
assess their contributions to the systematic uncertainty. 

6 Conclusions 

The branching ratio Br(b — ^ DDX) has been measured using an inclusive joint probability 
method with data collected by the OPAL detector at LEP. The result 

Br(b ^ DDX) = (10.0 ± 3.2(stat.)l2:9(det.)l9V(pliys.))% 

*The momenta of the D hadrons are scaled by the maximum possible momentum for D hadrons produced in 
continuum e"'"e~ annihilations at ^/s = 10.58 GeV. 
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is consistent with the inclusive measurement of Br(b — > DDX) by DELPHI: Br(b DDX) = 
(13.6 lb 3.0(stat.) ± 3.0(syst.))% [6]. Several significant sources of systematic uncertainty that 
are investigated in this analysis were not assigned uncertainties in the DELPHI analysis. These 
sources partially account for the difference in the size of the estimated systematic uncertainties. 

Combining this Br(b DDX) result with previous experimental determinations of Br(b — > 
charmonium) = (2.4 ± 0.3)% [15] and Br(b^ no charm) = (0.7 ± 2.1)% [6] in Equation 1 
gives He = 1.12lo;io- This value is consistent with the average value of ric = 1.166 ib 0.033 [2], 
from previous measurements carried out at the Z^. The measurements of ric and Br(b — >■ DDX) 
obtained in this analysis are consistent with theoretical calculations [1]. 
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